Exponential Estimates in Adiabatic Quantum Evolution^ 

Alain Jovef^ and Charles-Edouard Pfisterj^] 
Abstract 

We review recent results concerning the exponential behaviour of transition probabilities 
across a gap in the adiabatic limit of the time-dependent Schrodinger equation. They range from 
an exponential estimate in quite general situations to asymptotic Landau-Zener type formulae 
for finite dimensional systems, or systems reducible to this case. 

2P ' 1 Introduction 

Q\ ■ 

Q\ • The notion of adiabatic evolution or adiabatic process is an important theoretical 

concept, which occurs at several places in Physics. In Quantum Mechanics, this process 
is usually described by the equation ih-^jt/j(t') = H(et')tp(t'), where s is a small pa- 
rameter, such that 1/e gives the typical time-scale over which the Hamiltonian changes 
significantly. Setting h = 1 and introducing a rescaled time, t = et', we can rewrite the 
time-dependent Schrodinger equation for the evolution operator U as 



CO 



d 

ie—U(t,s)=H(t)U(t,s), U(s,s)=l, Vt, s G R. (1) 



The adiabatic limit corresponds to the singular limit e — > of the equation (|l|). If the 
O ■ state of the system is an eigenfunction ip(t ) for the eigenvalue e(t ) at t = t , then, 
in the adiabatic limit, the state of the system at time t — t\ is an eigenfunction ip{t{) 



for the eigenvalue e(ti), provided the energy-level e(t) is isolated in the spectrum of the 
Hamiltonian H(t) for all t in the time-interval [to>^i] [4]. When the system performs a 
cycle, H(t ) = H(ti), and the energy-level e(t) is non-degenerate, the eigenfunction ip(ti) 
defined in the adiabatic limit differs from ij)(to) by a phase, which can be decomposed 
into an ^-dependent dynamical phase and an ^-independent geometrical phase related to 
the spectral subspaces visited during the adiabatic evolution. This is the fundamental 
observation of Berry, which gave rise to extensive developments (see the collection of 
papers in [32]). In this note we review a complementary aspect, namely the estimation 
of the probability that the final state of the system is not an eigenstate ip(ti) for e(ti); 
such a transition is called adiabatic or nonadiabatic transition in the physical literature. 
There are two kinds of results. On the one hand for small systems (or systems reducible 
to this case), typically two-level systems, one can derive explicit formulae for the proba- 
bility of such a transition. A particular case is the Landau-Zener formula used when the 
two levels display an avoided crossing. For general systems, on the other hand, one can 
usually obtain upper estimates for the transition probability only. Those estimates are 
reviewed in the second part. The whole discussion is lead in a scattering setting under 
the assumption that the time- dependence of the Hamiltonian is analytic. More precisely 
HI The self-adjoint family {H(t)} t <=-R, is defined on a common dense domain D of a 
separable Hilbert space Ti; it is uniformly bounded from below; for each <fi G D, H(t)<p 
has an analytic extension in a fixed open strip ScC including the real axis. Moreover, 
there exist if (±oo), two self-adjoint operators defined on D, and constants C, a > such 
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thatfort>0, sup slt+iseS \\(H(t + is)-H(±oo)) ( p\\<cM+^^Pl. 

H2 The spectrum o~(t) of H(t) consists in two disjoint parts o~(t) = o~i(£) U o~2{t), such 
that inf tg R dist (cx(t), cr 2 (t)) > g > and cri(t) is bounded. 

Denoting by Pj(t) the spectral projector associated with a jit), j = 1,2, the transi- 
tion probability between the subspaces Pi(to)TC et -P^i)^ is defined as V(ti,to,e) = 
\\P2(t 1 )U(t 1 ,t Q )P 1 (t )\\ 2 . The Adiabatic Theorem of Quantum Mechanics [4], [21], [28], 
[1] states that V(ti,t ,e) = 0(e 2 ) as soon as H(t) is C 3 . Under the assumption HI 
(analyticity an d existence of H(±oo)) we have "P(+oo, — oo, £:) = 0(e°°) [25], [29], [22]. 
The difficult questions are to prove exponential estimates and asymptotic formulae for 
V(+oo, — oo, e). 

2 Asymptotic formulae 

Assume the Hamiltonian H (t) is an n x n hermitian matrix with non degenerate spec- 
trum a{t) = {ei(t), ■ ■ ■ , e n (t)}. We determine up to a global phase factor a corresponding 
basis of normalized eigenvectors (fj(t), j — 1, • • • , n, by requiring that (<pj(t)\<p'j(t)) = 0, 
teK. We can expand the solution ip{t) = U(t)ip(0) of ([]]) as 

m = T,^t)e-'Io^/^t), (2) 

3=1 

where the c'^s are complex valued coefficients and the phase factors e~ l Io e ^ t " ,dt ^ £ are 
introduced for convenience. It follows from HI that the limits lim t ^± 00 Cj(t) = c., (±oo) 
exist. Thus, choosing Cj(— oo) = 5jk as initial conditions, we get that |cj (oo) | 2 yields 
the transition probability from the eigenspace associated with ek(— oo) to the eigenspace 
associated with ej(+oo), which we denote by Vjk(s)- Born and Fock [4] showed that 

Cj(t) = Cj(-oo) +0(e) (3) 

uniformly in t G R, from which the Adiabatic Theorem Vjk{s) = 0(e 2 ) follows. 

In case n = 2, for real symmetric analytic Hamiltonians, the pioneering works [24], [34] 
and [5] gave arguments in favour of the Landau-Zener-Dykhne formula V^is) — e~ 27 / £ 
with 7 > explicited below. A convincing derivation of this formula based on the 
integration of ([!]) in the complex t-plane can be found in the important contribution [9]. 
However, for complex hermitian matrices, this formula misses a prefactor of geometrical 
nature, as we now show. The proof is based on the fact that the solution ip of (|I|) 
is analytic throughout the strip S, since H is, whereas <pj and Cj have multivalued 
extensions. More precisely, their only possible branching points are the complex crossing 
points zo G S, such that ei(zo) = e2(^o)- Generically, ei{z) — e\{z) ~ y/z — z Q . Let rj be 
a negatively oriented loop based at the origin, which encircles Zq. Denoting by f{z\r]Q) 
the analytic continuation along i] of a function f(z) defined in a neighbourhood of the 
origin, we can write ei(z\rj ) = e 2 (z). It follows that ipi(z\r}o) is proportional to ^2(2) 
and we define #2i(?7o) ^ C by (/^(z^o) = e^ 21 ^ (p 2 (z) . Since il>(z\r) ) = ip(z), we deduce 
from the foregoing and (§) the key identity 

Cl ( z | r?0 )=e^^)e l 4 ei(y)d2 ' /£ C2 (^). (4) 

Let us assume, for simplicity, that there exists a unique generic crossing point zq in 
S with Im zq > 0. By analyticity of ip, we can integrate ([!]) from —00 to +00 with 
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Cj(— oo) = Sij along the real axis or along any path (3 C S which passes above z . 
Denoting by C\ the analytic continuation of C\ along such a path /3, we have, due to 

(i); C2(+oo) = e^ 1021 ^ ^ l ho ei( - 2 ^ dz ^(^(+00). It thus remains to prove the estimate 
ci(+oo) = 1 + 0(e), which is the equivalent along (3 of the estimate (|]). As is well 
known from complex WKB methods, such an estimate can be proven under a global 
dissipativity condition on the path j3 only. By definition, a path (3 C S parametrized 
by t 1— > i3(t) with limtioo Re (3{t) = ±00 is called dissipative (for the indices {1,2}) if 
Im fp, t Jei(z) — e 2 (z))dz is non-decreasing in t G R, where means integration from 
— 00 to (3(t) along /3(s), —00 < s < t. Sumerizing the above discussion we have 
Theorem 1 [12] Assume HI, H2 for H(t) a 2 x 2 matrix and suppose there exists a 
unique generic crossing point z Q in S with Im z > 0. Then, provided there exists a 
dissipative path in S going from —00 to +00 above Zq, we have for e > small enough 

V 21 {e) = e 2Im <M*) e 2Im k ^ d ^[i + {e)). (5) 

Remarks: The complex factor e -16 * 21 ^ ^ is actually obtained by analytic continuation 
along i]q of the quantity yielding the geometric phase when considered on the real axis. 
When H(t) is real symmetric, Im ^21(^0) — an d we recover the Landau- Zener-Dykhne 
formula. However, in general, Im ^21(^0) 7^ 0, thus giving rise to a non-trivial prefactor. 
This prefactor was independently discovered by Berry [3]. 

The (non-trivial) questions of existence of dissipative paths and competition between 
several crossing points in S are analyzed in details in [12]. 

In case no dissipative path above Zq exists, an exponential bound follows from integration 
along dissipative paths close to the real axis [12], [20]. See also Theorem 2. 
We refer the reader to tables 1 and 2 in [16] and references therein for similar asymptotic 
formulae in more general 2-level systems as well as in non-generic situations. 

When H(t) is an n x n matrix with n > 2 the same strategy is applicable in principle. 
Let us see this for n = 3, to fix the ideas, in the following simple setting. Suppose there 
exist two distinct crossing points z and z 1 only in S, with Im Zj > 0. Let z be a generic 
crossing point for the analytic continuations of e\(t) and e2(t) around a loop 770 encircling 
zq only, as above. Assume the analytic continuation of e${t) in a neighborhood of 770 is 
analytic at zq, as is generically the case. By hypothesis, z\ is a generic crossing point for 
the analytic continuations of e2(t) and e^(t) around a similar loop r\\ encircling z\ only, 
whereas the analytic continuation of e\{t) in a neighbourhood of r\\ is analytic at z\. We 
get, as above, that (U) holds and must be completed by 

C2 (z\ Vl )=e ie ^e I ^ z ' )dz ' /s C3 (z), (6) 

where 832 {vi ) is defined similarly. Assume the composition of the two negatively oriented 
loops, based at the origin, rjo followed by 7/1 can be deformed into one negatively oriented 
loop, based at the origin, encircling both zq and Z\. Then, if c\ denotes the analytic 
continuation of c\ along a path going from —00 to +00 above zq and Zi, we get from 

(1) and © c 3 (+oo) = e-^ 1 ^)e" i ^o ei(2 ' )d2 ' /£ e- ie32 ^)e" i ^i e2(2 ' )d2 ' /£ ci(+oo). Again, it 
remains to prove ci(+oo) = 1 + 0(e) with initial conditions Cj(— 00) = 8\j. Such an 
estimate can be proven in a dissipative domain, defined as follows. A domain P C S is 
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dissipative if it extends from — oo to +00 and if any point z G P can be reached from 
—00 by a dissipative path @ 2 C P for the indices {1, 2} and by a dissipative path /3 3 C P 
for the indices {1,3}. However, generally, dissipative domains going above Zq and Z\ do 
not extend from —00 to +00 [6], [9]. The way out is to perform a perturbative analysis. 

Assume H(t, 5) depends on a supplementary parameter 5 > 0, satisfies HI for each 
fixed 5 and is regular enough in (t, 5). Suppose that for 5 = 0, the (analytic) eigenvalues 
ej(t,5 = 0), t G R, labelled in increasing order at t = —00, are non degenerate, except 
for two generic real crossings at t < t\ where ei(t , 0) = e 2 (to, 0) and ei(ti, 0) = e 3 (ti, 0). 
Assume these degeneracies are lifted as soon as 5 > 0, turning t and t± into avoided 
crossings. As a consequence, for 5 > small enough, there exist crossing points z (5) 
and Zi (<!>), with the properties discussed above and lim^ Zj(S) = tj. Furthermore, the 
existence of a suitable dissipative domain P can be proven perturbatively so that we get 
Theorem 1' [11] Assume H(t,5) is a 3 x 3 matrix satisfying the hypotheses (loosely) 
stated in the previous paragraph. Then, for 5, e > small enough, 

-p u ( £ ) = e 2Im0 21 (r, o ,<5) e 2Im J VQ ei{z' ,5)dz> / e ^Im e 32 ( m ,8) g 2Im e 2 (z>,6)dz> /e ^ + Qt £ \\ tj\ 

where 0(e) is uniform in 5 and lim^o ^ 8) = lim^o Im f Vk e j{ z ' 1 S)dz' = 0. 

Remarks: The above result and generalizations in [11], [18] show that it is possible to 
get asymptotic formulae for certain transition probabilities in n-level systems displaying 
avoided crossings, despite the difficulties induced by the notion of dissipativity. 
The same approach was used when two levels only display one avoided crossing [10]. It 
yields back Theorem 1 without need to check dissipativity conditions in the complex 

plane. Moreover, assuming the generic behaviour e 2 (t) — ei(t) ~° \/a 2 t 2 + 5 2 , we get the 
Landau-Zener formula: 2 Im J* e\{z' , 5)dz' /e = —^(1 + 0(5)) and Im #21(770, 8) = 0{5) 
[10]. Further refinements of the Landau-Zener formula can be found in [27]. 

3 Exponential estimates 
Let us come back to general abstract time-dependent Hamiltonians. In such cases we 
can prove that the transition probability across the gap is exponentially small. 
Theorem 2 [13] Let H(t) be a time- dependent Hamiltonian satisfying assumptions HI 
and H2. Then, there exist C, 7 > and Eq > such that for all < e < Eq, 

V(+oo,-oo,e) = i lim o ||P 2 (tiMti,to)Pi(to)ir < Ce^' £ (8) 

The proof is made by integrating (|l|) in the complex plane along generalized dissipative 
paths suitable for operators (see [13], section 5), following [12]. A weaker result can be 
found in [19]. Similar estimates were known for finite dimensional systems of ODE's, 
see e.g. [7] and the literature quoted in [6]. Theorem 2 was then recovered by different 
methods, some of which allowing a better control on the exponential decay rate 7 as a 
function of the gap g. See [26] for an approach using microlocal analysis, [30] and [15] 
for superadiabatic techniques, see below, and [33] for a pseudo-differential point of view. 
It follows from [15] and [26] that 7 > eg, g large, for some constant c. 

4 Superadiabatic renormalization and reduction theory 

The root of this method is the work [8] and subsequent generalizations and adapta- 
tions [29], [31], [14] on iterative schemes. Since the work [21], the Adiabatic Theorem is 
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often proven by showing that the adiabatic evolution V(t, s) defined by the equation 

ie^V(t,s) = (H(t)+ie[Pl(t),P 1 (t)])U(t,s), V(s,s)=l, Vt, s G R (9) 

satisfies the intertwining property V(t, s)Pj(s) = Pj(t)V(t, s) and \\U(t,s) — V(t,s)\\ = 
0(e) [1]. Such an approach can be used in a wider context, see [23]. An iterative 
scheme consists in generating from H(t), by standard perturbation in e, a sequence of 
Hamiltonians {H q (t, e)} gG N with the following features: 

i) The Hamiltonians H q (t,e) share the same general properties as H(t). In particular, 
the gap hypothesis is true for H q (t,e) so that the spectral projectors Pj >q (t,e) are well 
defined and tend to Pj(t) as e — > 0. 

ii) The adiabatic evolution V q (t,s), defined by (|9|) with H q (t,e) and Pj <q (t,e) in place of 
H(t) and Pj(t), approximates U(t,s) up to order 0(e q ) instead of 0(e). 

It follows that the transition probability between spectral subspaces of H q (t,e) is 
0(e 2q ) instead of 0(e 2 ) for finite times. Berry [2] formally showed for two-level systems 
that one can push the estimates up to exponential order by truncating the scheme at 
optimal q G N. In the present setting, Nenciu proved such exponential estimates in 
[30] using closely related ideas. The result of Nenciu was recovered in [15] by optimal 
truncation of the iterative scheme proposed in [14]. 

Theorem 3 [30], [15] Assume HI and H2. Then there exist constants C, 7,£o > and 
a self adjoint operator H*(t,e) defined on D such that for all < e < Eq and t G R, 
\\H*(t,e) — H(t)\\ < C*£:/(l + |i|)( 1+a ). Moreover, the adiabatic evolution V*(t, s) associated 
with H*(t,e) satisfies sup ts£K \\V*(t, s) — U(t, s)|| < Ce" 1 ^ as well as V*(t, s)Pj t *(s,e) = 
Pj t *(t,£)V*(t,s). Finally, 7 > Cg, for g large. 

Remarks: Theorem 2 becomes a corollary of Theorem 3, since P,>(t,e) -5 Pj{±oo). 
The evolution V*(t, s) is called superadiabatic evolution due to the exponential estimate. 

The main interest of this construction, however, is that it allows to set up a rigorous 
reduction theory. Assume Pi(t)7i is finite dimensional, say two-dimensional, such that 
a t(f) = { e i(t)i e 2(^)} where ej(t) satisfy the gap assumption. Then the transition proba- 
bility V2i{s) between eigenspaces associated with e\{— 00) and e2(+oo) can be computed 
modulo errors C(e~ 7/ ' £ ) by replacing the evolution U (t, s) by the superadiabatic evolution 
V*(t,s) which describes the evolution inside the Pj^(t,e)/H. This leads to an effective 
two-dimensional problem, with a corresponding effective 2x2 Hamiltonian shown to be 
a perturbation of H(t)Pi(t) and to have the same analyticity properties as H(t) in S. 
Hence, provided the conditions on the analytic continuations of Cj stated in Theorem 1 
are satisfied, we get with the same notations [15] 

V 21 {e) = e 21iae '^e 2lm ^o ei{z)dz/£ (l + 0{e)) + 0{e^ /£ ) (10) 

The condition 7 > |2Im / ei(z)dz\ is satisfied if g is large enough [15] or in the avoided 
crossing regime described in Theorem 1' [10], [11]. 

Another use of the superadiabatic renormalization consists in performing the analysis 
of section 2 for finite dimensional systems when the decomposition (||D is replaced by 

ijj(t) = YJj=i Cj,*(t)e~ l ^° ej '*^ ,£ ^ dt ^ 6( fj,*(t, e) where e^i, e) and (fj t *(t,e) are the eigenval- 
ues and eigenvectors of H*(t,e). After showing that they have analytic continuations in 
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S with similar properties as those of ej and ipj, we get very accurate asymptotic formulae 
by making use of Theorem 3 "in the complex plane". Under the conditions of Theorem 1 

V 21 (e) = e 2ime 21 „( £ ) e 2Im 4 ) e M(^ £ )^/ £ (1 + 0(e-^))), 7 > 0, (11) 

as shown in [17], [11]. The loop r] and 9 2 \^{e) are defined as in section 2. 
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